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Topological quantum computation is a promising technique to achieve large-scale, error-corrected 
computation. Quantum hardware is used to create a large, 3-dimensional lattice of entangled qubits while 
performing computation requires strategic measurement in accordance with a topological circuit 
specification. The specification is a geometric structure that defines encoded information and fault-tolerant 
operations. The compilation of a topological circuit is one important aspect of programming a quantum 
computer, another is the mapping of the topological circuit into the operations performed by the hardware. 
Each qubit has to be controlled, and measurement results are needed to propagate encoded quantum 
information from input to output. In this work, we introduce an algorithm for mapping an topological 
circuit to the operations needed by the physical hardware. We determine the control commands for each 
qubit in the computer and the relevant measurements that are needed to track information as it moves 
through the circuit. 

Quantum computing is at a pivotal point in development. The major enabling factor is the greatly 
enhanced understanding of highly efficient quantum error- correcting codes which allow reliable opera- 
tion of a quantum computer even when its individual qubits have very large error-rate. This paper 
focuses on the generic class of measurement-based topological quantum computers' (TQC), which combine 
several advantages that make them one of the most promising paradigms to build scalable quantum computers in 
the near future^ TQC is based on advanced topological quantum error- correcting codes, where the quantum 
hardware produces a well defined array of qubits that are arranged and connected using quantum entanglement 
links in a regular 3D lattice. 

TQC allows, in principle, a reliable realisation of large quantum algorithms, such as the Shor algorithm to factor 
numbers in polynomial time''' '^, a task that is not efficiently possible on a classical computer. However, imple- 
menting complex functionality requires specific design flows based on clear abstraction levels and automated 
synthesis steps. In case of classical circuits, the design stack typically includes algorithmic level, register-transfer 
level, gate level, transistor level, and mask level (sometimes, more levels are defined). In over 50 years of classical 
circuit design, a large variety of automated and semi-automated methods to translate a design from a higher level 
to a lower level have been developed. 

For quantum algorithms based on the TQC paradigm, a design stack can also be defined (Figure 1). A quantum 
algorithm, such as the Shor algorithm, is first transformed into a quantum circuit consisting of quantum gates 
(each gate is a manipulation of typically 1, 2 or 3 qubits that can be described by a unitary matrix). This step is 
required for TQC and all other quantum architectures, and substantial body of work on gate-level design and 
optimisation has been published both in Quantum Physics and in Computer Engineering communities"'"^''. 

TQC circuits allow for a generic geometric description that abstracts from details, in particular the underlying 
error correction, which is sufficient to capture relevant algorithmic interactions. The geometric structures in the 
TQC circuit, representing encoded qubits and encoded gates are embedded (via physical qubit measurements) 
within the physical 3-dimensional lattice. The model provides direct support for only a restricted subset of 
quantum operations, due to symmetries necessary for error correction to function, and if the high level circuit 
contains arbitrary gates, they must first be mapped to the supported operations using a series of well-known 
steps'^'^'-^"*. After that, the circuit can be translated into its canonical geometric representation where each gate has 
a corresponding region of the 3D lattice where it is implemented. The first automated optimisation method for a 
restricted case of two-dimensional descriptions was presented in^^. 

This paper focuses on the next synthesis step: mapping of the geometric description to the physical level. The 
output of the mapping consists of two parts: A direct instruction set for the quantum hardware and a classical 
instruction set for data processing. The geometric description of the TQC circuit specifies a layout of how logical 
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Figure 1 | The TCQC design stack. At the lattice and Hardware levels, 
there are two elements to the design stack; one classical and one quantum. 
The quantum component is the instruction sets that tells each qubit in the 
computer how it is measured. The classical component is keeping track of 
measured data in order to successfully propagate information from input 
to output in the topological circuit. 

qubits should be arranged within the lattice of physical qubits. The 
physical size and separation of these structures is proportional to the 
error correction strength needed for computation. TQC operates by 
performing measurements on physical qubits within the lattice. The 
type of measurement depends on whether a given physical qubits is 
inside or outside the geometric region defining a logical qubit. The 
instruction set for the quantum hardware consists of translating the 
geometric structure to a specific set of individual measurement 
instructions for the physical qubits in the lattice. 

The classical data processing information is also required for 
byproduct calculation. Within a TQC circuit, large sets of physical 
qubits form what is known as a correlation surface. These sets con- 
nect the inputs of a quantum circuit to the outputs. The overall parity 
of these correlation surfaces (when they are measured) tells us how 
the logical information must be corrected in order to successfully 
propagate from input to output. The classical control software of a 
TQC must know these structures in order to successfully operate. 
Therefore, for a TQC circuit specification, the determination of the 
physical qubit sets that define correlation surfaces and how they 
interact is non-trivial and constitutes the by far most complex part 
of the mapping algorithm. 

The result of the mapping algorithm is a 3D lattice of qubits along 
with instructions how to use it. This is vaguely related to the general- 
purpose hardware and the customised software of a classical 
microprocessor. Note that TQC is technology-independent: any 
implementation where qubits can be arranged in two or three dimen- 
sions can be used to realise this model*"'" '^. Therefore, our work closes 
the gap between the high-level synthesis (which reaches up the geo- 
metric description level) and technology- oriented physical imple- 
mentation approaches. It is the foundation for the first complete 
automated flow from the quantum algorithm to the implementation. 

The physical qubits are arranged in a large three-dimensional 
lattice and are entangled to form a single, massive, universal 
quantum state. A universal state is a specific quantum state that 
can be used to implement any fault- tolerant quantum algorithm 
via measurement of each qubit in the state^'. The unit cell of the 
lattice is illustrated in Fig. 2a). Prior to actual computation, each 
qubit is initialised in the superposition state, | +) = (|0) -|- 1 1)) 
/ \/2. Qnce all qubits are initialised, a two-qubit gate is used to 
entangle any two qubits that are connected in the lattice with an 
edge. The unit cell is the building block of the lattice, and its structure 
is repeated along three dimensions, which are denoted by w (width), 
h (height) and t (time). TQC is implemented via sequential measure- 





a) 

Figure 2 | Structure of the Topological lattice, (a) Complete 3D lattice 
cell. For the cell (at lattice coordinates cc e CEL) qubits marked black are 
Sea 2nd gray qubits are F^^ (b) A primal cell (marked gray) and eight dual 
cells (physical qubits and entanglement are not represented). 

ment of a 2D cross section of the full 3D lattice. Information is 
propagated along the third axis of the lattice via the entanglement 
links, with processing occurring during these teleportation steps. 
Hence the third dimension of the physical lattice is identified with 
the temporal axis (t-axis). 

One general unit cell of the lattice has 27 (3 X 3 X 3) vertex 
positions, 18 of which are occupied by physical qubits, denoted by 
gray and black circles in Figure 2a). A lattice with mc„ X mch X mCf 
unit cells has a total of (2mc„ + 2) • (2mCf, + 2) • (2mCf +2) positions, 
because neighboring cells share a side. All physical qubits in the 
complete lattice form the set TQCC. The positions not occupied by 
qubits form the set CEL. A position (cvf, ct, ch) e CEL must be the 
center of a unit cell. The physical qubits at the faces of such a cell (six 
gray circles in Figure 2a) are denoted by Fe„ ,.h,cn ^nd the qubits at the 
sides (12 black circles in Figure 2a) are denoted by S^w.ch.ct- The set of 
all qubits of a cell is Ccw,ch.ct '■ = Ecw,ch,ct U Scw,ch,ct ■ Note that TQCC is 
the union of these sets for all the cells in the lattice. 

The topological lattice contains two self-similar lattices that are 
interlaced (see Figure 2b)). The primal lattice (with corresponding 
C, F", y) contains all the cells having odd coordinates, and the dual 
lattice (with corresponding C', F*, S'') contains all the cells having 
even coordinates. If eight dual cells of the lattice are arranged into a 
cube, then a primal cell exists at the intersection point, and therefore 
the set of all physical qubit coordinates in the lattice can be expressed 
as TQCC = FP{JF''. 

Physical qubits (organised in cells arranged in a 3D-lattice, as 
described above) are used to encode logical qubits. The underlying 
lattice is a computational resource for the logical layer, similar to how 
the hardware is a computational resource for the physical layer. 
However, at the logical layer, all operations are error-corrected, 
implying that logical qubits are initialised, operated on and measured 
using logical operations. It should be noted that the logical informa- 
tion is not decoded prior to any logical operation and afterwards re- 
encoded; operations are applied directly on encoded physical qubit 
states. Valid operations within the logical space are designed specif- 
ically such that physical errors do not cascade to cause uncorrectable 
faults. Therefore, the TQC model is known as a fault-tolerant model 
of quantum computation. 

The geometric structures, which define encoded qubits, are cre- 
ated by measuring relevant physical qubits in the Z basis (|0, 1) 
states). This effectively creates a hole in the lattice, known as a defect. 
It is the defects that encode information. Two types of logical qubits 
can be defined depending in which of the lattices the Z measurements 
of physical qubits are performed and consequently where defects are 
created. The physical qubits sets C f and D"* C P' represent the 
physical qubits that are internal to defects which are measured in the 
Z basis, with D''nS'' = 0andD''nS'' = 0.Primaldefectsarecreated 
by measuring qubits in the set D*' and Dual defects are created by 
measuring qubits in the set D''. 

Figures 3a) and b) illustrates the structure of two defects which 
encode a single logical qubit (the unit cells that are not relevant for 
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Figure 3 | Qubit sets used for information processing. Figure (a) is for loop operators (The X (Z)-eigenstate information for primal (dual) 
defects) and Figure (b) is for chain operators (TheZ(JO-eigenstate information for primal (dual) defects). These sets define defects (green circles), logical 
operators (orange squares) and correlation surfaces (blue triangles) for a single logical qubit undergoing the identity operation. A shorthand notation is 
illustrated in the inserts. This figure illustrates for minimum sized defects (and hence the smallest distance quantum code). Figure (c) illustrates 
initialisation in the Z, Xbasis and state injection. State injection is where in individual qubit is measured in a rotated Zbasis (the vertex of the pyramids) 
and then expanded to form an encoded defect, this allows us to achieve universality. The superimposed black vertices and lines are used to illustrate graph 
edges and vertices in the mapping algorithm. 



the defects are not shown). In this structure a temporal axis is defined 
running from the front of the structure to the rear and this encoded 
qubit is undergoing a simple identity operation, i.e. the input of this 
topological circuit starts on the 2D layer at the front of the image and 
via measurement is teleported to the 2D layer at the rear. There are 
three sets of physical qubits that are relevant for the mapping. One is 
the sets D''*'" C P"''* that are used to define primal (dual) defects. 
These qubits are illustrated in Figures 3a) and b) as green circles and 



are measured in the |0, 1) basis to simply remove them from the lattice. 
The second set are known as logic operators. Logic operators are sets of 
physical qubits whose parity determine the logic state of the encoded 
information. Two types of logic operators exist, one is a loop operator 
that completely encircle an encoded defect. The second is a chain that 
connects both defects (illustrated by orange squares in Figure 3b). 

For primal (dual) defects, the loop operators correspond to the 
logical X (Z) operator while the chain corresponds to the logical Z (X) 
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Figure 4 | Evolution of correlation surfaces during a braided CNOT gate, (a) The geometry of a CNOT gate where the dark (dual) qubit acts as a control, 
(b) X'' correlation surface on the control (black lines) is intersected by the target defect (green line) creating a correlation surface, Xi" in accordance with the 
properties of the CNOT. (c) illustrated that the Z' surface of the control does not become correlated, as required. 



operator. The final set of relevant qubits are illustrated as blue trian- 
gles in Figures 3a) and b). These sets are known as correlation sur- 
faces. The correlation surface connects logic operators from an input 
layer in the cluster to an output layer. The ring and chain operators 
need to be connected from the input layer (front of the image) to the 
output (rear of the image). The qubits illustrated in blue perform this 
operation. If each qubit in blue is measured in the | + , ^) state (X 
basis) we connect the chain (ring) operators from the input side to 
the chain (ring) operators on the output side. If we wished to measure 
the logical state, by measuring and calculating the parity of the logical 
operator we would measure each physical qubit in the Z basis. By 
measuring these qubits in the X basis instead, the logical information 
is teleported to the next layer in the lattice, with a correction that is 
determined by the parity of the X basis measurements. The correla- 
tion surfaces form either sheets or tubes. For primal (dual) defects, 
sheets are the logical Z'' {X') surfaces and tubes are the logical X'' (Z^) 
surfaces. In order to track information as it moves through the lattice, 
we need to know which qubits are associated with these surfaces, as 
the parity of these measurements will dictate the correction that 
needs to be applied to encoded information as it moves through 
the lattice. 

The primary goal of this work is, given a geometric definition of a 
TQC circuit, to determine the information needed by the (classical) 
control software of a quantum computer in order to operate TQC 
hardware that implements a circuit. This information includes the 
measurement basis of the physical qubits. In Figures 3a) and b), the 
solid circle qubits (defect-internal) are measured in Z basis whereas 
the triangle qubits (correlation surfaces), square qubits (input/out- 
put) and white qubits (qubits for error correction) are measured in X 
basis. Moreover, the measurement outcomes for qubits in correlation 
surfaces are needed to calculate corrections to encoded data as they 
propagate through more complicated topological structures that are 
not considered here. Therefore, the mapping procedure must 
identify the physical qubits for each correlation surface (X*"' C P" 
and Z^-'' C 

Other geometric structures, corresponding to specific logic opera- 
tions, are shown in Figure 3c). For a primal pair of defects, the 
horseshoe structure corresponds to initialising a qubit in the 1 0) state, 
the pair of parallel defects correspond to initialisation in the | + ) state 
and the third structure corresponds to state injection. State injection 
allows us to introduce an arbitrary state into the TQC lattice by 
measuring a single physical qubit that exists at the vertices of the 
pyramid structure. State injection is instrumental in the TQC model 
to achieve universal computation^"^^. In Figure 3c) the temporal axis 
runs from left to right. In order to perform logical measurement, we 
use the exactly same structures, but in the reverse temporal direction 
(this allows for measurement in the |0, 1) and | + , ^) basis, i.e. Z and 
X bases). The regions that are marked as grey and green correspond 
to the sets of physical qubits whose parity defines the parity of the 
initial encoded qubits. 

The concept of state injection in the topological model is instru- 
mental for achieving universality. The power of quantum computa- 



tion is intricately related to the introduction of these injected states 
and their consumption to realise non-Clifford gates through tele- 
portation. However, in the context of this mapping algorithm, there 
is essentially no difference in the classical computation associated 
with the geometric structures in Figure 3c). While the right most 
structure in Figure 3c) represents the injection of a non-trivial 
encoded state, the sets of physical qubits that define byproduct 
operators for this newly created state are no more complex to cal- 
culate than the sets associated with standard Clifford initialisation. 
Therefore, the introduction of injected logical states and the realisa- 
tion of non-Clifford gates in a topological circuit poses no complexity 
problem for the mapping algorithm. 

The mapping problem of a TQC circuit is the accurate determina- 
tion of physical qubits in the topological lattice corresponding to 
relevant correlation surfaces that define the operation of the circuit. 
The input to this algorithm is a suitably optimised topological circuit 
that is compatible with the error-corrected model and hence the 
underlying computational hardware. Algorithmic constraints at 
both the physical and logical level have already been taken into 
account via the construction and optimisation of the input circuit. 
Before formally introducing the problem, we illustrate it through the 
example of a cnot operation implemented via a braid operation 
between two encoded defects (Figure 4). A correlation surface that 
begins on primal defects cannot terminate on dual defects and vice 
versa. Therefore, as defects are braided, a correlation surface extends 
around the qubit involved in the braid. It is this interaction that 
entangles the logical qubits. 

Figure 4 illustrates how correlation surfaces are defined, using the 
shorthand notation that is illustrated in the inserts of Figures 3a) and 
b). In Figure 4b) we show the X'' correlation surface (sheet) on the 
control qubit (a dual type defect). Due to a braid with a primal defect, 
this correlation surface must extrude around the second defect and 
forms a tube. This correlates the outputs of the two qubits when the 
control qubit is in an X-basis eigenstate. Figure 4c) shows the Z'' 
correlation surface (tube) on the control qubit. Here, the surface is 
not perturbed by the braid. Therefore the output remains uncorre- 
cted when the control is in a 2-basis eigenstate. This circuit corre- 
lates output qubits in precisely the way needed to realise a CNQT 
operation'*'^^. 

These surfaces are not specified with the TQC circuit, and there- 
fore the mapping procedure must derive them from the geometric 
structure and map them to the actual set of qubits that are used in the 
lattice such that they can be tracked as the circuit is implemented. 

Results 

Problem formalism. The mapping procedure takes a geometric 
description of a TQC computation as input and calculates a data 
structure Q which contains the required information of each logical 
qubit. Each logical qubit in Q is specified by a set, ct, of directed 
segments, given by their end points {begin, end) e CEL^. Recall 
that CEL is the set of lattice coordinates that are not occupied by 
qubits and can be centres of unit cells. Each segment has one of six 
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Figure 5 | Example of the algorithm illustrated in Eqn. 4. The above 
illustrates for the graph = ({A, /}, {{A, B), (B, C), ... (J, A)}). For 
example, starting from the vertex start = A, the first possible operation is 
reduce{E, F) and a first sub-sheet is found SUBS = {(£, G)}. Furthermore, 
because C, D, G, Hare co-linear will result in remove(D) and remove(G). 
The cycle is traversed until the start is reached again, and the first traversal 
completes. After the second traversal neither reduce nor remove were 
applied. Therefore, the reshape(A, B, C) is applied, and a second sub-sheet 
is found SUBS = {{A, C), (£, G)}. Finally the last two sub-sheets are 
inferred leading to SUBS = {(A, i), (C, I), (A, C), (£, G)}. The complete 
sheet is found by combining aU the sub-sheets according to the algorithm 
in Eqn. 3. 

possible directions {±w, ±h, ±t}.To convert t/ to Q, an intermediate 
representation y where each logical qubit is described by a directed, 
cyclical graph, is calculated first and then used to derive Q. 

Geometric description and representation hy graphs. For each logical 
qubit o", a directed cyclical graph Gg = {K^, V^,) is defined. The set K 
of possible graph vertices contains representations of unit cells from 
set CEL. There is a bijective mapping coord between K and CEL: 
coord{k) yields the three (w/h/t) lattice coordinates of the cell cor- 
responding to k.Kc^K represents the end points of the segments in 
ff, and Vg represents the connections between the end points. A 
graphical example for this construction is offered in Figure 3c). 

Mapping algorithm. Representation of Logical Qubits by Set Tuples: 
The tuple representation = [1,D'^jI,01jI,xI,Z'^) of a logical 
qubit a contains all the information which is needed to implement 
it and its interactions with other qubits on a quantum computer. The 
ultimate outcome of the mapping procedure described next is the set 
Q which contains tuples for all logical qubits a. 

The specifics of the mapping algorithm can be found in the 
Methods section. The edges in V„ form a Hamiltonian path, contain- 
ing all vertices from K^. During the mapping, the path can be either 
modified by remove(a); fl e K that removes vertex and its incident 
edges, and creates a new edge between its neighbours, or the insert{a, 
b, c); a,b, c e K function that deletes the (a, c) edge, and creates the 
new edges (a, b) and {b, c). 

Within q^, I stands for the type of the logical qubit (primal or dual). 
D'^ includes all defect- internal physical qubits (i.e., those to be mea- 
sured in the Z basis, whereas all other qubits will be measured in the X 
basis). ll and O'^ are physical qubits that define inputs and outputs, 
respectively. In Fig. 3, the qubits in D'^ were shown in green and the 
qubits in and C/^ in orange. Sets and include all physical 
qubits that are part of the X and Z correlation surface, respectively. 
Finally, jI is the set of injection points. Note that the primal (dual) 
qubits are defined using coordinates on primal (dual) cells. 

As mentioned earlier, the Z correlation surface is a sheet for a 
primal logical qubit and a tube for a dual logical qubit. The X cor- 
relation surface is a tube for a primal logical qubit and a sheet for a 
dual logical qubit. For convenience, we introduce functions sheet: 
Q^TQCC which returns Z^ for primal logical qubits and for 
dual logical qubits and tube: Q^TQCC which returns the corres- 
ponding tube in a similar way. 



Mapping examples. The presented algorithms were implemented 
and their execution will be illustrated by examples. The example 
from Figure 5 illustrates how a complete sheet is computed. A 
more complex example, for the computation of tubes and sheets, is 
presented in Figure 6, and procedure starts from the TQC 
computation described by the geometry presented in subflgure a). 
The circuit consists of 3 logical CNOTs performed between 4 logical 
qubits: 3 primals (the qubits 1, 2, 4) and a dual (qubit 4). The circuit 
will implement the following sequence: CN0T{3, 2), CN0T{3, 4), 
CN0T{3, 1). 

The quantum circuit in Figure 6a) is converted to a geometric form 
appropriate for the TQC model (Figure 6b)). The mapping will first 
determine the tube surfaces of each logical qubit and the results 
shown in Figure 6c) will be obtained. Figure 6d) contains the result 
of the physical qubit sets forming sheet surfaces. For the same dual 
logical qubit (the third qubit in Figure 6a)) all the intermediate steps 
of the sheet computation are graphically presented in Figure 6e). 

Algorithmic complexity. The algorithmic complexity of the map- 
ping procedure is analysed in the following. For the computation of 
all the qubit sets, computed by the algorithm in Eqn. 1, the mapping 
procedure requires a single qubit-cycle traversal, while the compu- 
tation of the coordinates is straightforward (see the illustration 
presented in Figure 3a and the corresponding legend). Therefore, 
for a given logical qubit q„, the runtime complexity is linear in the 
number of vertices (\K„\) of the graph-cycle representing the 
geometry. 

For the runtime complexity analysis of the sheet mapping proced- 
ure a worst-case geometry has to be defined. Such a geometry, when 
mapped in the 3D, will have to require a maximum number of 
reshape-mle applications in order to be able to compute the corres- 
ponding sub-sheets. The search of a geometry wiU consider that the 
algorithm from Eqn. 4 randomly selects a vertex from the cycle- 
graph (Line 1), and that the vertex is used as a pivot for the 
reshape-operation. These assumptions imply, that even after a 
sequence of reshape operations is applied, there is no possibility to 
find a sub-sheet and to reduce the number of vertices, and thus a 
further reshape is necessary. For a graph with vertices, a worst- 
case situation arises when \Kg\ — 3 vertices are arranged in a pattern 
similar to the one in Figure 7, where the red cube indicates the start 
vertex and the cycle is traversed clockwise. The graph vertices (where 
edges are joined) are not represented. 

After a first traversal of the graph, the reduce operation could not 
have been applied, and a reshape followed by corresponding remove 
operations will transform the graph like in the right panel of Figure 7. 
It can be noticed, that until the step-like construction is not reshaped, 
no reduce can be applied. However, after each reshape, co-linear 
vertices can be removed, thus minimising \Kj\. Overall, the complex- 
ity of the sheet-finding procedure is bounded by ©(jK^rj^). 

While the mapping algorithm itself scales polynomiaUy with the 
size of the cycle graph, we cannot make any claims regarding the size 
of the cycle graph given an optimised topological circuit. The size of 
the cycle graph is essentially related to the number of 90° "pivots" 
made by defects in the topological circuit we do not know yet how 
big, in general, this set is for a given circuit specification. Current 
methods in systematic optimisation for topological circuits are still in 
their infancy^" and determining the scaling of the number of these 
pivots for a given number of logical qubits is unknown. However, the 
physical size of the topological circuit (in terms of the physical 
volume of cluster needed to implement the circuit at a fixed error 
correction strength) has marginal influence on the scaling of the 
mapping algorithm itself Instead, the scaling of the mapping algo- 
rithm wiU be dominated by the total number of pivots made by the 
input circuit specification. 

Correctness. Checking the correctness of the algorithms, implies 
verifying their termination and the fact that the correct physical 
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Figure 6 | Defining tubes and sheets, (a) A standard circuit identity to perform a CNOT between two primal defects', mediated by a single dual defect that 
acts as control for the three braiding operations, (b) The geometric description of the circuit employed. The gray areas indicate the braiding of the logical 
qubits; (c) The sets of qubits forming tube correlation surfaces for the topological structure: red qubits form the set of the dual surfaces, while the primal 
tube surfaces are formed by the blue qubits; (d) The sets of qubits forming sheet correlation surfaces: the surface of the dual logical qubit is represented as 
the set of red physical qubits, while the sheets for the primal logical qubits are only illustrated, (e) Intermediate steps in finding the sheet correlation 
surface for the dual logical qubit with the algorithm detailed in the methods section. The edges of the graph used for computing the surface are indicated 
with thick lines. This figure is slightly rotated from Figure (c). 
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Figure 7 | Worse-case geometric structure for algorithmic efficiency. The 

above is an example of defect structures that require a maximum number 
of reshape-rule apphcations to calculate the corresponding sub-sheets. 

qubit coordinates are computed. Again, for the algorithm in Eqn. 1, 
the termination is guaranteed because a single traversal of the cycle- 
graph is necessary. The correctness of the computed coordinates is 
shown by comparing the algorithm instructions with Figure 3. The 
direction d (Line 4) is associated with the green line in the figure, 
where b and e are the CEL coordinates of the edge vertices. For 
example, let us consider that the two cells from the lower defect 
are b and e. This implies, that by selecting the two neighbouring 
coordinates (Line 7), the co-ordinates of the green marked qubits 
(the defect set D) is computed. The remaining 4 qubits (Line 8) that 
do not belong to D are the light blue marked qubits, which are 
associated with the tube correlation surface that surrounds the 
defect region. One has to note, that for defects with a cross-section 
larger than a cluster-ceU, the set union operations on Line 10 and 
Line 12 are to be interpreted as (A\{a}) U(({fl}\j4)n{fl})i meaning 
that if element a existed in the set A it wUl be removed, otherwise it 
will be included. 

The termination of the sheet-finding algorithm (Eqn. 4) can be 
illustrated by starting from the fact that the geometric description is 
mapped into a 3D representation where only 6 segment directions 
are possible (see the discussion of the reshape and reduce rules). A 
geometrically described defect configuration of a logical qubit (in the 



a) 




absence of any possibility to apply reduce or remove) will have an 
even number of edges in its associated graph. Considering all the 
worst-case geometries (an example is Figure 7), it can be seen that 
after a graph-traversal, either reshape or reduce (followed by remove) 
can be be applied. The number of maximum consecutive reshapes is 
bounded by \K„\ . Thus, the number of graph-vertices is continuously 
reduced, and the termination of the algorithm is guaranteed. 

The correctness of the coordinates computed from the sub-sheets 
is shown by comparing the Lines 3 and 4 of Eqn. with Figure 3b). For 
a logical qubit of type Z, the set of all possible side qubits (depicted 
gray in Figure 2 for a single cell) is S'. For the two defects from 
Figure 3b) the instruction on Line 3 will return the coordinates of 
both the orange and the blue marked qubits, which intersected with S' 
wiU return only the blue qubits. These are the ones necessary for the 
computation of a sub-sheet. If two sub-sheets overlap (which is pos- 
sible, given the way the algorithm in Eqn. functions), then the inter- 
section set is not part of the set of qubits defining the complete sheet 
of a logical qubit (Line 4). 

A final detail that needs explanation is that the approach, by find- 
ing sub-sheets, is correct. Initially, one has to consider that for TQC, 
by deforming the defect geometries, does not change the computa- 
tion, because the topology of the description is relevant, and not its 
geometry. For this reason, for each individual logical qubit described 
as a graph, the reshape operation is valid, because the sheets of two 
topologically equivalent geometries will differ only in sub-sheet. For 
example, in the last panel of Figure 8b), the subsheet BD is the 
difference between the equivalent geometries: before and after apply- 
ing the reshape operation. 

Another fact about the correctness of sub-sheet-finding appears 
when one considers the example presented in Figure 8a). Let us 
consider, the three qubits A, B, C e S'. Deforming the defect implies, 
that the qubits U, V are not part of the defect anymore, and that D, E 
e S' are parts of the new boundary. This is equivalent to stating that 
the sub-sheet containing the qubits D, E was created, and that after 
the change the complete sheet will contain ({A,B,C} U{£',£} a S') . 




b) 




Figure 8 | Using defect geometries to define correlation surfaces, (a) Changing the geometry of a defect. The trajectory of the defect (indicated by the 
dashed line) modifies the boundary of the sheet from the initial boundary A, B, C, extending the boundary A, B, C involving the qubits U, V, Wand the 
resulting boundary A, D, E. (b) Three types of graph operations to define a surface. The first, applying reduce( C, D) results in \KJ being decreased by 2, 
because C, Dare removed. After the operation, because A, B, E, Fcorrespondto co-linear lattice coordinates, remove{B); remove{E) can be further applied, 
and \K„\ again decreased by 2. In the second case, the reduce(C, D) is applied. The mirrored vertices V"'' = [C, D'] are computed, with coord(D') = 
coord(B), thus V"' = {C, B\. The vertices from VZ'' = {C',B}\{B,£} = {C'} wiU be inserted, and the vertices from V]^, = {B} will be deleted. In the last 
case, the effect of the reshape(B, C, D) operation is that vertex C is replaced by vertex C. Afterwards reduce{A, B) and reduce(D, E) can be applied. 
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It should be noted that the correctness of the mapping algorithm 
already assumes that the given topological circuit has been verified 
with respect to the initial circuit specification. In general, a com- 
pressed topological circuit may bare little resemblance to the original 
un-optimised circuit. Before a topological circuit is input to the map- 
ping algorithm, we assume that it has been verified against the initial 
specification and consequently the qubits sets generated from the 
mapping algorithm will faithfully realise the original circuit when 
implemented in hardware. Optimisation of topological circuits may 
employ a technique known as bridge compression"" which connects 
closed loop topological structures together. The algorithm detailed in 
this work is not currently designed to handle circuits containing 
bridged components, however modifying the algorithm to process 
arbitrary structures will be possible as arbitrary geometric shapes can 
still be decomposed into simple cycle-graphs. These modifications 
are the subject of future work. 

Discussion 

In this paper we presented the first mapping algorithm that translates 
a topological quantum circuit into the specific measurement patterns 
required by the quantum hardware. We have also illustrated the 
necessary algorithm for the construction of correlation surface 
needed to appropriately track information as it propagates through 
the circuit. The mapping of a topological circuit specification to the 
physical measurements in the lattice is one of many classical compo- 
nents necessary for realisation of quantum circuits. This classical 
library consists of both office and online components. Offline soft- 
ware processing consists of the compilation and optimisation of a 
topological quantum circuit'^'^" ™ and the mapping of this structure 
to the physical lattice. Online classical processing relates to the real- 
time decoding and processing of error correction data using algo- 
rithms such as minimum weight perfect matching^" and renormali- 
sation group techniques''. Ultimately, all classical processing will 
need to be integrated together to form a complete software package 
for controlling topological quantum computers. 

This work is not only a necessity for the future operation of topo- 
logically protected quantum computation, but it is also a require- 
ment for appropriate verification packages for topological circuit 
synthesis. The accurate construction of correlation surfaces allow 
us to compare the operation of a topological circuit against the 
abstract quantum circuit it is derived from. This is important as 
optimisation techniques for topological quantum circuits often result 
in structures that bare little resemblance to the original circuit spe- 
cification. The further development of verification protocols is the 
focus of future work. 

Methods 

The mapping procedure. The algorithms to compute the classical control 
information are presented in detail in this section. For this purpose, we first describe 
several functions to translate information between the geometric description, the 
cycle-graphs and the tuple representation of logical qubits. 

Function type: V ^ {init, measure, inject, defect] returns the type of a geometric 
segment that was translated into a graph edge. The segment type is used after physical 
qubits associated with the segment have been calculated in order to decide which set 
from q„ they belong to. 

Function set takes q„ and the segment type as inputs and returns the corresponding 
coordinate set of physical qubits that belongs to segments oiq^ with the specified type. 
For example, physical qubits from a segment of type inject \i will be added to the set f^, 
and for init segments the is returned. 

Function dir. V ^ [±w, ±h, ±t} yields the lattice-direction of the segment 
represented by the given edge. A further function ngh"{k), k ^ K, neZ is used to 
indicate the n-th neighbour of vertex k in the direction of the Hamiltonian path. 

Function coord: CEL returns the vector representing the lattice coordinates of 
a graph vertex k ^ K. The inverse function coord^^ is used to infer the vertex given the 
lattice coordinates. 

The function mirr. K; mirr{a) — coord^^ (coord{ngh^^{a)) + coord{ngh{a)) — 
coord(a)) returns vertex a mirrored at the line through its predecessor and successor 
in the cycle-graph, while the cist: —> K; clst(a, b, c) returns either b or c depending 
which is closer to vertex a. 



The mapping starts by constructing for each geometrical description of a logical 
qubit (7 from Q the cycle-graph G^g3^. Then, the corresponding tuple q^^Q is initi- 
alised by setting / to primal or dual according to the desired type and 
D',=ll = Oi=}l=Xl = Zi = 0). 

Sets i^, 7^, 0|^, 7^, tube{q„) are constructed by the algorithm shown in Eqn. 1. 
Recall that tubeiqc^) is either or Z^, depending on whether q„ is primal {/ — p) or 
dual (/ — d). The algorithm traverses the cycle-graph G^j. The lattice- direction d 
associated with each edge of the graph is computed. For each cell cc of a segment, F^c is 
its complete set of physical face qubits (see Section). Two face qubits along the 
segment are defect- internal and added to D^c- The remaining four qubits are part of 
the tube correlation surface and added to set r„ that is then added to tube{q„). 
Coordinates of injection points are always found at the middle of the {b, e) segment 
(Line 15). The coordinates of the qubits on iniilmeasure segments are added to the 
corresponding sets l„ and as given by set{q„, type) at Line 12; 

Require : q^ - (/, D, 7, O, ]. X, Z)\q„& Q\ ae Q 

Require :G^-(7:^,y^),G^G 3^ 

1 : stori'*— random ke 

2: ck<— start 

3 : repeat 

4: d'^dir{{ck, ngh{ck))) 

5: b*— coord (ck); e<—coord{ngh{ck)); 

6: for all cc G (/?, e) along d 

7: Dcc*-{p\p — cc± 1 along d, p G Fee} 

8: r,,<-P,,\7)cr 

9: if type{{ck, ngh{ck))) ^ defect (1) 
10 : tube{qfj)<- tube{q„) U T^c 

11: end if 

12: set{qa, type)*— set (q^, type)\JDcc 
13: end for 

14: if type{{ck, ngh{ck))) — inject 

15: }^}yj{{b-\-e)/2} 

16: end if 

17: ck<—ngh{ck) 

18: until start ^ck 

19: return q^ 

The sheet surfaces of each G^^ are found by a procedure that iteratively reduces the 
graph until it consists of just two vertices. The procedure constructs a sheet by finding 
one sub-sheet in each iteration, where a sub-sheet contains physical qubits bounded 
by a rectangle in either wh, wt or ht plane. Each sub-sheet is specified by two points 
{ss\ ss^); 55' e CEL; where 55' define the diagonal of the rectangle and have one equal 
component. Sub-sheets found by the procedure in different steps are disjoint, and the 
union of all found sub-sheets is the complete sheet. Figure 6d) illustrates a complete 
sheet composed of sub-sheets, and Figure 6e) shows the progressive calculation of 
sub-sheets (both figures will be discussed in more detail later). 

The procedure transforms the graph by eliminating or moving vertices, while sub- 
sheets are calculated. Without giving a full formal proof for correctness of these 
transformations (such a proof will require advanced quantum computing concepts 
which are out of this paper's scope), we illustrate the modification of a sheet boundary 
in Figure 8. The dashed line in the figure is indicating the direction of one defect 
involved in generating the sheet surface SHEET. By changing the direction of the 
defect, the boundary A, B, C e SHEET is transformed into A, D, E e SHEET'; 
SHEET = SHEET U {D, £}. In terms of sub-sheets, the sub-sheet {B, E) is added after 
the change of defect direction. 

The algorithm illustrated in Eqn. 4, used to compute the SUBS set of sub-sheets for 
each logical qubit, wUl remove vertices of the graph, by constantly traversing it (Lines 
3, 20), until only 2 vertices are left (Line 2). Eqn. 3 takes each SUBS set and constructs 
sheet{(T) of each logical qubit tuple o. For the subsheet {ss^, ss^) the function 
boundingboXss = x A, (} [fnin(^ss]^ ,max(^ss^^^ r\TQCC constructs the set of 
coordinates able to represent physical qubits. 

Considering the structure of the lattice and the way the geometrical descriptions 
are constructed, it can be noticed that the cycle-graphs have an even number of 
vertices. It should be also noted that the insert operation is not directly applied by the 
algorithm, and the number of vertices is modified only during the reduce operation, or 
after the reshape operations is applied. 

The reduce operation. The algorithm illustrated in Eqn. 4 modifies the graph by 
applying the the reduce operation. The path is reduced if for 3 consecutive edges, the 
first and the second edge represent opposite associated directions into the lattice (Line 
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10). For example, in Figure 8b) this is the case for the edges {B, C), (C, D), (D, E), where 
(B, C) and (D, E) have opposite directions. 

The reduce{a, b); a, b ^ K operation is defined as the sequential application of: 

Rm — {a,b}] 

na ^ngh-^ (a) ; - ngh^ (b) ; Ng - { J 
V'^'^ — {dst{a, mirr{b), na),clst{h, mirr{a), n\,)} 



(2) 



removeiy) ; for all vGRm 
insert{na, Vj, n\,)\iox aU VjGy^^^; 
remove{v^)\ for aU e ^^^f ' 

We illustrate the reduce operation by applying it to vertices C and D of Figure 8b) . The 
sets Rm - {C, D}, JV^ - {B, £}, - {d5f(C, mirr{D), B), clst{D, mirr(C), £)} are 
constructed. Because mirr{C) = E and mirr{D) = B, the set V^""* ~ {B, E} is equal to Ng 



yre 



- 0. After the vertices from Rm are removed, no further vertices are 



inserted or removed because the corresponding sets are empty. However, for the 



example in Figure 8b) this is not the case as | V^^c 
removing the vertex B and inserting vertex C. 



= 1, thus effectively 



The reshape operation. In order for the graph to be reduced, it may be required to 
represent an equivalent geometrical description. Thus, vertices are not removed or 
deleted, but moved (jiCo^j remains constant). The re5/jiipe((a,fo,c); a,/?, c e Koperation 
is the sequential application of: 

remove{h)\ insert{a, mirr{b), c)-. 

In the context of the algorithm in Eqn. 4, the function is called if during a complete 
graph traversal the reduce operations could not be applied. For the example of 
Figure 8b), where reshape{B, C, D) is called, the resulting graph will be obtained by 
removing C and inserting C ~ mirr{C). Applying reshape for a second time at the 
same position would undo the initial application, thus reduce{B, C, D) is the inverse 
of reduce{B, C, D). Therefore, the start pivot (Line 3), used for checking if a traversal 
completed (Line 20), is updated (Line 24). 

The maximum number of cycle traversals is 0{\KJ^), when after each traversal a 
reshape operation is required, while the number of consecutive reshape operations is 
bounded by |K] — 3. 

The number of vertices of the cycle-graph is initially even due to the shape prop- 
erties of the geometrical descriptions. Moreover, it remains even during the execution 
of the algorithm. Operation reduce eliminates exactly two vertices from the graph (set 
Rm), while the sets of further added and deleted vertices {V^^ and Vj^f , respectively) 
are always of the same size. Operation reshape does not add or delete vertices. 
However, three consecutive vertices may represent a straight line after a reshape 
operation, in which case they are replaced by two vertices (Line 14). It can be shown 
that one further vertex elimination must follow, keeping the overall number of ver- 
tices even. 

Require -.SUBS 

Require -.q^ - (/, D, /, O, X, 2); q^GQ 
1: sheet{qtj)*- 
2: for all ssgUBS do 

(3) 

3: SSQ<— [q\coord{q) e boundingbox{ss)DS'} 
4: sheet{q,)^{sheet{q^)\JSSQ)\{{sheet{q,)r\SSQ)) 
5 : end for 
6: return a 



Require - (K,V) {otggQ 

1 : SUBS^ 

2: while |K| >2 do 

3 : start -e- random keK 

4; ck<— start 

5; compact <- false 

6: a*-ngh^{ck): b*-ngh^{ck)) 

7 : repeat 

8; if dir{{cKa))^ -dir((b,ngh(b))) then 

9; compact true 

10: reduce{a,b) 



SUBS^SUBSU{ngh{ck),a){J{ngh{ck),b) 
else 

if dir{{ck,a)) — +dir{{a,b)) then 
remove{a) 
compact*— true 
else 

ck<-a 
end if 
end if 
until stark — ck 
if compact —false then 

reshape (start ,ngh{start) ,ngh^ {start)'j 
SUBS^SUBSU {start, ngh^ (start)) 
start*— ngh(start) 
end if 
end while 
return SUBS 



(4) 
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